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Abstract. New facilities to measure neutrino-nucleus cross sections, such as those 

■ possible in conjunction with spallation neutron sources, could provide an experimental 
C*~) \ foundation for the many neutrino-nucleus weak interaction rates needed in supernova 

models. This would enable more realistic supernova models and provide a greatly 
f~>) . improved ability to understand the physics fundamental to supernovae by comparison 

t-H ' of these models with detailed observations. Charged- and neutral-current neutrino 

■ interactions on nuclei in the stellar core play a central role in supernova dynamics 
and nucleosynthesis as well as being important for supernova neutrino detection. 

' Measurements of these reactions on judiciously chosen targets would provide an 

^p" 1 . invaluable test of the complex theoretical models used to compute the large number 

of neutrino-nucleus cross sections that are needed. 

. 

i PACS numbers: 97.60.Bw, 25.50,+x, 23.40.-s 

• i-H , 

ft 

1. Core Collapse Supernovae 

Core collapse supernovae are among the most energetic explosions in the Universe, 
releasing 10 46 Joules of energy in the form of neutrinos of all flavors at the staggering 
rate of 10 57 neutrinos per second and 10 45 Watts. Marking the death of a massive star 
(mass > 8 — 10 solar masses) and the birth of a neutron star or black hole, core collapse 
supernovae are a nexus of nuclear physics, particle physics, fluid dynamics, radiation 
transport, and general relativity. They serve as laboratories for physics beyond the 
Standard Model and for matter at extremes of density, temperature, and neutronization 
that cannot be produced in terrestrial laboratories. The 10 44 Joules of kinetic energy 
and rich mix of recently synthesized elements delivered into the interstellar medium by 
the ejecta of each supernova make core collapse supernovae a key link in our chain of 
origins from the Big Bang to the formation of life on Earth. 

The center of a massive star as it nears its demise is composed of iron, nickel, and 
similar elements, the end products of stellar nucleosynthesis. Above this iron core lie 
concentric layers of successively lighter elements, recapitulating the sequence of nuclear 
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burning that occurred in the core during the star's lifetime. Unlike prior burning stages, 
where the ash of one stage became the fuel for its successor, no additional nuclear energy 
can be released by further fusion of the maximally bound, iron peak nuclei. No longer can 
nuclear energy production stave off the inexorable force of gravity. When the iron core 
grows too massive to be supported by electron degeneracy pressure, the core collapses. 
This collapse continues until the core reaches densities similar to those of the nucleons 
in a nucleus, whereupon the repulsive core of the nuclear interaction renders the core 
incompressible, halting the collapse. Collision of the supersonically falling overlying 
layers with this stiffened core produces the bounce shock, which drives these layers 
outward. However, this bounce shock is sapped of energy by the escape of neutrinos 
and nuclear dissociation and stalls before it can drive off the envelope of the star (see, 
e.g., P). The failure of this prompt supernova mechanism sets the stage for a delayed 
mechanism, wherein the intense neutrino flux, which is carrying off the binding energy 
of the proto-neutron star, heats matter above the neutrinospheres and reenergizes the 
shock 00] • The heating is mediated primarily by the absorption of electron neutrinos 
and antineutrinos on the dissociation-liberated free nucleons behind the shock. This 
process depends critically on the neutrino luminosities, spectra, and isotropy, i.e., on 
the multigroup (multi-neutrino energy) and multi-angle neutrino transport between the 
proto-neutron star and the shock. 

Realistic supernova models will require extremely accurate neutrino radiation 
hydrodynamics, but, as we will show, unless advancements in neutrino transport and 
hydrodynamics are matched by equally important advancements in nuclear and weak- 
interaction physics, the accuracy of the former must be called into question. For electron 
and electron-neutrino capture rates, and in general all neutrino-nucleus cross sections, 
detailed shell model computations must replace the parameterized approximations that 
have been used in the past. For example, recent work |3] on electron capture up to 
atomic mass 65 has shown that these parameterized rates can be orders of magnitude in 
error. The electron capture rate is dominated by Gamow- Teller resonance transitions, 
with the Gamow- Teller strength distributed over many states. Previous parameterized 
treatments [3] place the resonance at a single energy, and this energy is often grossly 
over- or underestimated relative to the Gamow- Teller centroid computed in realistic 
shell model computations, leading in turn to rates that are often too small or too 
large, respectively. (If the centroid is underestimated, more electrons can participate 
in capture. The reverse is true if the centroid is too high.) To take full advantage of 
these improved rates, improvements in the tracking of the nuclear composition will be 
needed. The use of a single representative nucleus, which has been the standard until 
now in virtually all supernova models [HI El, can underestimate the rates of electron and 
electron- neutrino capture during the critical core collapse phase (see, e.g., jH]). 

While generations of nuclear structure models will afford ever greater realism 
in the calculation of stellar core properties and the interactions of nuclei with the 
neutrinos flowing through the core, nuclear experiments must be designed and carried 
out that will serve as guide posts for the theoretical predictions that are required to 
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produce the large number of rates that enter into any realistic supernova or supernova 
nucleosynthesis model. In particular, we must have neutrino-nucleus cross section 
measurements that will help gauge the accuracy of neutrino capture and scattering 
predictions important during stellar core collapse and during a-rich freezeout and the 
r-, p- and v nucleosynthesis processes after core bounce. 

2. Spherically Symmetric Supernova Models 

Although four decades of supernova modeling have established the textbook explanation 
discussed in the previous section, fundamental questions about the explosion mechanism 
remain. Is the neutrino heating sufficient, or are multidimensional effects such as 
convection and rotation necessary? Can the basic supernova observable, an explosion, be 
reproduced by detailed spherically symmetric models, or are multidimensional models 
required? Many observations of core collapse supernovae show clear violations of 
spherically symmetry. For example, neutron star kicks jU] and the polarization of 
supernova emitted light ^Hj cannot arise in spherical symmetry. Nonetheless, it remains 
uncertain whether such observations point to fundamental features of the supernova 
mechanism or merely side effects. To answer this question, simulations in one, two, and 
three dimensions must be coordinated. 

The neutrino energy deposition behind the shock depends sensitively on the 
neutrino luminosities, spectra, and angular distributions in the postshock region. 
Ten percent variations in any of these quantities can make the difference between 
explosion and failure in supernova models jTH H2j. Thus, accurate multigroup 
Boltzmann neutrino transport must be considered in supernova models. Past spherically 
symmetric simulations have implemented increasingly sophisticated approximations 
to Boltzmann transport: simple leakage schemes ^3], two-fluid models [Hj, and 
multigroup flux-limited diffusion [T5J [THJ H3 EH]- A generic feature of this last, most 
sophisticated approximation is that it underestimates the isotropy of the neutrino 
angular distributions in the heating region and, thus, the heating rate fHJ|2n]- With 
these limited transport approximations came the possibility that the failure to produce 
explosions in the past may have resulted from the incomplete neutrino transport. 

To address this possibility, the ORNL theoretical astrophysics group has used 
Boltzmann neutrino transport to successfully simulate the core collapse, bounce and 
evolution for more than 500 milliseconds after bounce of a number of spherically 
symmetric models. Even with our implementation of complete Boltzmann neutrino 
transport, via AGILE-BOLTZTRAN, models in both the Newtonian [21] and the general 
relativistic j22j limits have failed to produce explosions for a range of progenitor masses 
from 13-40 solar masses (SSlEHESj- These results are supported by the failure of a 15 
solar mass model to explode using a different implementation of Boltzmann transport 
[2E]- In all cases, neutrino heating propels the shock to a radius of more than 150 km; 
however, the heating is insufficient to eject the envelope, so the shock recedes. In the 
general relativistic model starting from a 40 M & progenitor |24j . the formation of a black 
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Figure 1. We trace the shock fronts for Newtonian (left) and General Relativistic 
(right) supernova models which fail to produce explosions for a range of progenitors 

[23 El ESI 

hole is evident approximately one half second after bounce. 

Thus, we are beginning to answer some fundamental questions in supernova theory. 
In the past, it was not clear whether failure or success in supernova models was the 
result of inadequate transport approximations. However, because our failed explosions 
were obtained using accurate spherically symmetric hydrodynamics j^Zj, state of the 
art neutrino transport |28j . industry standard neutrino-matter interactions 16 and the 
equation of state of [7], our conclusion is unambiguous. Accurate neutrino transport 
alone does not overcome the failure of supernova simulations that assume spherical 
symmetry to produce explosions. This would suggest that changes in our input physics 
(i.e., weak interaction physics and nuclear physics in the Equation of State) and/or 
initial conditions (i.e., stellar evolution models) are needed and/or that multidimensional 
effects such as convection, rotation, and magnetic fields are required ingredients in the 
recipe for explosion. 

3. Neutrino-Nucleus Interactions and Core Collapse Dynamics 

While neutrino interactions with shock heated nucleons are the major source of the 
neutrino heating which drives the delayed shock, neutrino interactions with nuclei are 
also important before and after the passage of the shock. During the collapse of 
the stellar core, electrons are captured on heavy nuclei and free protons in the core, 
producing electron neutrinos that initially escape, deleptonizing the core. Where the 
shock forms in the stellar core at bounce and how much energy it initially has are 
set by the amount of deleptonization in the core during collapse. Deleptonization 
would be complete if electron capture continued without competition, but at densities 
of order 10 11-12 g cm -3 , the electron neutrinos become "trapped" in the core, and 
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Figure 2. The energy scales and composition as a function of density in the collapsed 
stellar core at bounce for a 15 M Q progenitor |2*5) . 

the inverse reactions — charged- current electron-neutrino capture on iron-peak nuclei 
and free neutrons — begin to compete with electron capture until the reactions are in 
weak equilibrium and the net deleptonization of the core ceases on the core collapse 
time scale. The equilibration of electron neutrinos with the stellar core occurs at 
densities between 10 12-13 gem -3 . Figure El summarizes the thermodynamic conditions 
throughout the core at bounce and displays the temperature, electron fraction (Y e ), 
electron chemical potential (// e )> and mean electron neutrino energy (E Ue ) in MeV as 
functions of the matter density. Also shown is the representative nuclear mass (A). 
The kinks near 3 x 10 7 gcm -3 mark the transition to the silicon shell. As the stellar 
core densities increase, the characteristic nuclei in the core increase in mass, owing to 
a competition between Coulomb contributions to the nuclear free energy and nuclear 
surface tension. For densities of order 10 13 g cm -3 , nuclei with mass > 100 dominate. For 
densities exceeding 10 14 gcm -3 , heavy nuclei are replaced by nuclear matter. Figure H3 
demonstrates that the nuclear composition shows a wide spread in mass, with species 
with significant concentrations having masses that differ by 40. Further, Figure H3 also 
shows that the abundances of nuclei with mass greater than 100 are significant as early 
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Figure 3. Details of the composition at two points during stellar core collapse 



as 10 11 gcm~ 3 . Thus, cross sections for charged- current electron and electron-neutrino 
capture on nuclei at least through mass 120 are needed to accurately simulate core 
deleptonization and to accurately determine the postbounce initial conditions. 

Improved shell model calculations of weak interaction rates for electron capture, 
positron capture, and (3~ and (3 + decays on nuclei relevant for stellar evolution 
(45 < A < 65) have become available in recent years EUJ. Heger et al have 
utilized these new weak reaction rates to improve upon the stellar evolution simulations 
of Woosley & Weaver (WW95)[24j, replacing the weak interaction rates for electron 
and positron captures and (3~ and (3 + decays. The WW95 models used the electron 
capture rates of Fuller, Fowler, & Newman (FFN) j3j and older sets of beta decay rates 

\\V2\. The most noticeable effect of these improvements is a marked increase in the 
electron fraction (Y e ) throughout the iron core before collapse. Because the final size of 
the homologous core, and therefore the shock formation radius, is proportional to the 
square of the trapped lepton fraction (YJ 2 ) at core bounce jSH], the persistence of these 
initial differences in Y e throughout collapse could have a discernible effect on the shock 
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energetics. 

In an initial attempt to determine the influence of these improved stellar evolution 
models on iron core collapse, Messer et al [31] performed parameterized radiation 
hydrodynamics collapse simulations using AGLE-BOLTZTRAN j23 EH1 HUH EH] • These 
simulations show no difference in initial shock formation position between the two sets 
of progenitor models when "standard" weak interaction physics is used during core 
collapse. Unlike stellar evolution simulations, the only available information about the 
nuclear composition in supernova simulations is that provided by the nuclear equation 
of state [7|: mass fractions of free nucleons, a-particles, and heavy nuclei, as well as the 
mean charge (Z) and atomic mass (A) of the heavy nuclei. This limited information 
precludes the use of reaction rates calculated for individual isotopes, which were, in 
any event, unavailable for A > 65. Figure El depicts the limits of the FFN and LMP 
data sets, demonstrating their usefulness for stellar evolution, but their inadequacy 
for core collapse. Instead, the Messer et al simulations used the standard neutrino 
emissivity from nuclei developed by Bruenn [TH] (see also [HE])- This prescription, based 
on the independant particle model, treats electron capture on heavy nuclei through a 
generic O/7/2 — > Of 5/2 Gamow- Teller transition in the average heavy nucleus identified 
by the equation of state. Because this treatment does not include additional Gamow- 
Teller transitions, forbidden transitions, or thermal unblocking [33 EDI E], electron 
capture on heavy nuclei ceases when the neutron number of the average nucleus exceeds 
40. As a result, electron capture on protons dominates the later phases of collapse. 
The steep dependence of the free proton fraction on changes in the electron fraction 
assures convergence of initially distinct Y e profiles to a similar Y e profile inside the 
homologous core and, consequently, to a similar shock formation radius. However, 
when the parameter in the models that controls the rate of electron capture on nuclei is 
adjusted so that the capture on protons does not dominate, significant changes in the 
core collapse dynamics are evident. This demonstrates the need for accurate treatment 
of nuclear electron capture during stellar core collapse. 

Recently, Langanke et al [12] have produced electron capture rates for a sample 
of nuclei with A = 66 — 112 using hybrid, shell-model-RPA calculations. As a major 
advance over the simple treatment of nuclear electron capture in previous supernova 
simulations, Hix et al [43 j have developed a treatment of nuclear electron capture based 
on these hybrid rates and the shell model electron capture rates from Langanke et al 
[1] (LMP) for A < 65. For the distribution of emitted neutrinos, the approximation 
described by Langanke et al [H] was used. To calculate the needed abundances of 
the heavy nuclei, a Saha-like NSE is used [H3 EE], including Coulomb corrections to 
the nuclear binding energy EH] , but neglecting the effects of degenerate nucleons 
|49j . This NSE treatment has been used in prior investigations of electron capture 
in thermonuclear supernovae |5()j . The combined set of LMP and hybrid model rates 
are used to calculate an average neutrino emissivity per nucleus. The full neutrino 
emissivity is then the product of this average with the number density of heavy nuclei 
calculated by the equation of state. With the limited coverage of rates for A > 65, 
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Enclosed Mass 

Figure 4. The electron fraction, entropy, density and velocity as functions of the 
enclosed mass at the beginning of bounce for a 15 M Q model. The thin line is a 
simulation using the Bruenn parameterization while the thick line is for a simulation 
using the LMP and hybrid reaction rate sets. 

this approach provided the most reasonable estimate of what the total electron capture 
would be if rates for all nuclei were available. This averaging approach also makes the 
rate of electron capture consistent with the composition returned by the EOS, while 
minimizing the impact of the limitations of the NSE treatment. 

This improved treatment of nuclear electron capture has two competing effects. In 
lower density regions, where the average nucleus is well below the N = 40 cutoff of 
electron capture on heavy nuclei, the Bruenn parameterization results in more electron 
capture than the LMP+hybrid case. This is similar to the reduction in the amount 
of electron capture seen in stellar evolution models [25J and thermonuclear supernova 
[oTI] models when the FFN rates are replaced by shell model calculations. In denser 
regions, the continuation of electron capture on heavy nuclei alongside electron capture 
on protons results in more electron capture in the LMP+hybrid case. The results of 
these competing effects can be seen in the upper pane of Figure which shows the 
distributions of Y e throughout the core at bounce. In addition to the marked reduction 
(11%) in the electron fraction in the interior of the PNS, this improved treatment also 
results in a nearly 20% reduction in the mass of the homologous core, which manifests 
itself as a reduction in the mass interior to the formation of the shock from .67 M & to .57 
M & in the LMP+hybrid case. A shift of this size is very significant dynamically because 
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the dissociation of .1 M & of heavy nuclei by the shock costs 10 51 erg, the equivalent of 
the explosion energy. There is also an 11% reduction in the central density and a 7% 
reduction in the central entropy at bounce, as well as a 10% smaller velocity difference 
across the shock and quite different lepton and entropy gradients throughout the core. 
In the outer regions, the higher electron fraction slows collapse, resulting in, for example, 
reductions of a factor of 5 in density and 40% in velocity in the vicinity of 0.8 M & . 

These differences in the behavior of collapsing stellar cores illustrates the 
importance of weak interactions with nuclei. At the onset of collapse, the nuclei of 
interest are clustered in mass between 50 and 70 along the neutron-rich edge of stability. 
Throughout collapse, decreasing electron fraction and increasing density pushes the 
composition to heavier and more neutron-rich nuclei, including nuclei 4-6 decays away 
from stability and with masses greater than 100. The KARMEN collaboration pioneered 
work in this regime, measuring the cross section for 56 Fe(z/ e , e~) 56 Co [5T] . which is one 
of the nuclei of interest early in collapse. However, the sheer number of potentially 
important species, and their instabilitymakes direct measurements of the rates an 
impossibility. Nonetheless, measurements of the relevant neutrino-nucleus interactions 
remain extremely valuable by providing the most relevant constraints on the theoretical 
models. The technique used by the KARMEN collaboration can be used to measure 
the electron-neutrino capture cross section on any of a wide range of nuclei, wherein the 
natural abundance is dominated by a single isotope and the element is a solid at room 
temperature. Several such nuclei are in the critical nuclear mass range: 55 Mn, 59 Co, 
89 Y, 93 Nb, 103 Rh and 115 In. Priority among these choices should be decided by their 
ability to constrain the theoretical rate calculations. 

4. Breaking Spherical Symmetry 

Multi-dimensional simulations allow investigation of the role convection, rotation, and 
magnetic fields may play in the explosion. Supernova fluid instabilities fall into two 
categories: (1) instabilities near or below the neutrinospheres, which we refer to as proto- 
neutron star instabilities and (2) convection between the gain radius and the shock, 
which we refer to as neutrino-driven convection. Proto-neutron star instabilities may aid 
the explosion mechanism by transporting hot, lepton-rich matter to the neutrinospheres, 
thereby boosting the luminosities at the neutrinosphere. Neutrino-driven convection 
may aid the explosion mechanism by boosting the shock radius and the neutrino 
heating efficiency, thereby facilitating shock revival. While a presentation of these 
macroscopic phenomena may seem out of place in a discussion of the role of microscopic 
neutrino-nucleus interactions, as we will demonstrate, the linkage of the microscopic and 
macroscopic scales are part of what make supernovae both interesting and challenging 
to study. 
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4-1. Proto- Neutron Star Instabilities 

Within the proto- neutron star (PNS), whose "surface" is defined by the neutrinospheres, 
a number of fluid instabilities may arise as the result of the lepton fraction and entropy 
gradients present. The lepton gradients are established by the deleptonization of 
the proto-neutron star via electron-neutrino escape near the electron neutrinosphere, 
whereas the entropy gradients result from the weakening supernova shock. (As the shock 
weakens, it causes a smaller entropy jump in the material flowing through it.) It is by its 
influence on these gradients that the microscopic physics links to the macroscopic fluid 
behavior. As an example, we return to our comparison in the previous section of models 
using an improved description of nuclear electron capture. Figure El shows the entropy 
and lepton gradients near the onset of PNS convection. The curves labeled BIO and 
LMPH10 show the results of the Bruenn and LMP+hybrid cases, both 10 ms after their 
respective maximum densities are reached. Clearly visible in the region between 20 km 
and 60 km are an inward displacement and steepening (flattening) of the negative lepton 
(entropy) gradient that results from our more accurate treatment of electron capture. 
The curves labeled LMPH16 present the LMP+hybrid model at 16 ms after bounce, 
when its shock has reached a radius matching that of BIO. Clearly the differences in the 
gradients described above are not transient, nor do they arise from the slower progress 
of the shock in the LMP+hybrid case. 

To illustrate the effects these changes in the entropy and lepton gradients could 
have on convection in the proto-neutron star, we have plotted the Ledoux criterion 
(d In p/d In Y e )(d In Y e /dr) + (d In p/d In s)(d In s/alr) in the lower pane of Figure El 
Instability in the Ledoux approximation arises when this expression is positive. The 
broad spikes on the right that coincide with the shock are due to the numerical 
spreading of the shock and, thus, are an artifact. In the model employing the Bruenn 
treatment, there is an extended region of instability between 30 and 60 km, while in the 
LMP+hybrid case, the unstable region is only 10 km across and displaced inward. This 
suggests that, with an accurate treatment of nuclear electron capture, PNS convection 
would be less extensive than previously thought and would occur deeper within the 
stellar core, though the similarity of the values of the Ledoux criterion suggests that it 
would occur with similar strength. 

While this spherically symmetric example demonstrates the importance of neutrino 
interactions in determining the dynamic behavior of the fluid, it is far from the 
final answer. The development of convection in the proto-neutron star is a radiation 
hydrodynamic phenomenon, rather than a purely hydrodynamic phenomenon. In this 
region of the stellar core, neutrinos and the stellar core fluid are coupled. The neutrinos 
have the ability to equilibrate an otherwise buoyant fluid element with its surroundings 
in both lepton number and entropy, rendering the fluid element nonbuoyant Thus, 
a determination of the degree of Ledoux convection in this region, and its extent, must 
await three-dimensional radiation hydrodynamics simulations. Past studies, in one 
respect or another, could not fully determine either. In one case [22], neutrino transport 
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Figure 5. The electron fraction, entropy and Ledoux criterion as functions of the 
radius soon after bounce. The thin line (BIO) shows results a simulation using the 
Bruenn parameterization at 10 milliseconds after bounce while the thick lines arc from 
the LMP+hybrid simulation set at 10 (solid, LMPH10) and 16 (dashed, LMPH16) 
milliseconds after bounce. 



was assumed to be spherically symmetric, which overestimates the degree to which the 
fluid element and its surroundings may equilibrate, and excluded the innermost 25 km of 
the core, in which it is unlikely that equilibration will occur owing to the large neutrino 
interaction cross sections and inefficient transport. On the other hand, other detailed 
studies |5H|IE%j did include this innermost region, but used radial-ray neutrino transport, 
which underestimates the degree to which equilibration will occur by disallowing lateral 
transport between rays. In the first study, vigorous convection did not occur, whereas 
in the latter, vigorous convection did occur and was confined to the inner 20 km of 
the core. The full vigor and extent will be determined by fully two-dimensional, and 
ultimately three-dimensional, models. 

Neutron fingers are another instability that may occur in regions of crossed gradients 
in lepton fraction and entropy. Neutron fingers are "doubly diffusive" instabilities, 
stemming from competing efficiencies of lepton number and entropy transport in the 
core. If, for example, a region of high entropy and low lepton fraction exists at a larger 
radius in a central gravitational field than a region of low entropy and high lepton 
fraction, and if entropy transport is more efficient than lepton transport, a fluid element 
at the interface of the two regions perturbed inward will sink, generating a low-entropy, 
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low-lepton-fraction (neutron-rich) finger. An early attempt to investigate these effects 
was made in the simulations by Wilson & Mayle ^H]- These one-dimensional models 
lifted spherical symmetry in an approximate fashion using a phenomenological mixing- 
length description for neutron finger convection inside the proto-neutron star, which 
boosted the neutrino luminosities, causing explosions. However, the strength of the 
neutron finger instability assumed by Wilson & Mayle is controversial, as shown by 
detailed numerical neutrino equilibration experiments |2J. 

This range of outcomes clearly demonstrate that to determine whether or not 
proto-neutron star instabilities exist and, if they exist, are vigorous will require 
simulations coupling three-dimensional, multigroup neutrino transport and three- 
dimensional hydrodynamics. Moreover, accurate neutrino interaction rates, for both 
the possibly unstable phase and the preceding collapse phase, will be essential. 

4-2. Neutrino -Driven Convection 

Neutrino-driven convection occurs in the region behind the stalled shock but above 
the gain radius as a result of the entropy gradient that forms as material infalls while 
being continually heated from below. Models show high-entropy, rising plumes and 
lower-entropy, denser, finger-like downflows beneath the shock. The shock is ultimately 
distorted by this convective activity. In the Herant et al [56J simulations, this large-scale 
convection led to increased neutrino energy deposition, the accumulation of mass and 
energy in the gain region, and a thermodynamic engine that ensured explosion. Herant et 
al used two-dimensional "gray" (neutrino-energy-integrated, as opposed to multigroup) 
flux-limited diffusion in neutrino-thick regions and a neutrino lightbulb approximation 
in neutrino-thin regions. In a lightbulb approximation, the neutrino luminosities and 
rms energies are assumed constant with radius. Herant et al described this transport 
scheme as the "most wanting aspect" of their simulations, stressing the need for 
more sophisticated multidimensional, multigroup transport in future models. In the 
Burrows et al |HZj simulations, neutrino-driven convection in some models significantly 
boosted the shock radius and led to explosions. However, they stressed that success 
or failure in producing explosions was ultimately determined by the values chosen for 
the neutrino spectral parameters in their gray ray-by-ray (one-dimensional) neutrino 
diffusion scheme. (In spherical symmetry (ID), all rays are the same. In a ray-by-ray 
scheme in axisymmetry (2D), not all rays are the same, although the transport along 
each ray is a ID problem. In the latter case, lateral transport between rays is ignored.) 
Focusing on the neutrino luminosities, Janka and Miiller [H], using a central adjustable 
neutrino lightbulb, conducted a parameter survey and concluded that neutrino-driven 
convection aids explosion only in a narrow luminosity window (±10%), below which the 
luminosities are too low to power explosions and above which neutrino-driven convection 
is not necessary. 

In more recent simulations carried out by Swesty 58 j using two-dimensional 
gray flux-limited diffusion in both neutrino-thick and neutrino-thin regions, it was 
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demonstrated that the simulation outcome varied dramatically as the matter-neutrino 
"decoupling point," which in turn sets the neutrino spectra in the heating region, was 
varied within reasonable limits. (The fundamental difficulty with gray transport schemes 
is that the neutrino spectra, which are needed for the heating rate, are not computed. 
The spectra are specified by a neutrino "temperature," normally chosen to be the 
matter temperature at decoupling. In a multigroup scheme, the spectra are computed 
differentially.) In two-dimensional models by Mezzacappa et al jSH]; the angle-averaged 
shock radii do not differ significantly from the shock trajectories in their one-dimensional 
counterparts, and no explosions are obtained. Neither the luminosities nor the neutrino 
spectra are free parameters. These two-dimensional simulations implemented spherically 
symmetric (ID) multigroup flux-limited diffusion neutrino transport, compromising 
transport spatial dimensionality to implement multigroup transport and a seamless 
transition between neutrino-thick and neutrino-thin regions. Recently, Fryer and Warren 
jHO] (using methods similar to [HE]) have demonstrated that three dimensional models 
exhibit convective behavior similar to the two dimensional models. This somewhat 
surprising preference for large scale convection in both two and three dimensional 
simulations is possibly explained by hydrodynamic instabilities in the stalled accretion 
shock j61j . Most recently, a simulation by [53], coupling 2D hydrodynamics with the 
ray-by-ray neutrino transport (performing independent calculations of the radiation 
transport along each radial direction) failed to produce an explosion. 

In light of the neutrino transport approximations made in all multidimensional 
supernova simulations to date, next-generation simulations will have to reexplore 
neutrino-driven convection in the context of three-dimensional hydrodynamic 
simulations that also implement more realistic multigroup three-dimensional neutrino 
transport. Furthermore, such simulations will only be as accurate as the neutrino matter 
interactions they include. In the previous section, we demonstrated the significant 
reduction in the density and velocity of the material infalling toward the shock due to 
an improved prescription for nuclear electron capture (see Fig. HJ). 111 addition, charged 
current neutrino capture (and neutral current inelastic neutrino scattering [02]) on heavy 
nuclei can alter the entropy and neutronization of this infalling matter prior to its arrival 
at the shock. Such changes in the pre-shock matter affect the shock propagation and 
thermodynamic conditions in the post-shock convective region, illustrating again the 
dependence of the macroscopic fluid dynamics on accurate microscopic physics. 

5. Supernova Nucleosynthesis 

Supernova nucleosynthesis is commonly divided into several "processes" , each of which 
is impacted by neutrino-nucleus interactions. (1) Explosive nucleosynthesis occurs as a 
result of compressional heating by the supernova shock wave as it passes through the 
stellar layers. In the inner layers of the ejecta, where iron group nuclei result from a-rich 
freezeout, interactions with neutrinos alter the neutronization, changing the ultimate 
composition. (2) Neutrino nucleosynthesis or the V process occurs due to neutrino- 
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induced nuclear transmutations in the outer stellar layers followed by shock heating. (3) 
The rapid neutron capture or "r" process may occur in the neutrino-driven wind that 
emanates from the proto-neutron star after the explosion is initiated. The neutrinos 
both drive the wind and interact with the nuclei in it. Early phases of this wind have 
also been suggested as the source of light p-process nuclei [S3]. Thus, neutrino- nucleus 
interactions are important to all core collapse supernova nucleosynthesis processes. 

5.1. Neutrinos and the a-Rich Freezeout 

One common property exhibited by recent spherically symmetric Boltzmann simulations 
PH is a decrease in the neutronization (which is equivalent to an increase in Y e ) 
of the inner layers of the ejecta due to neutrino interactions. This is a feature that 
current parameterized nucleosynthesis models can not replicate because they ignore the 
neutrino interactions. The neutronization of the ejecta is important because galactic 
chemical calculations and the relative neutron-poverty of terrestrial iron and neighboring 
elements strongly limits the amount of neutronized material that may be ejected into the 
interstellar medium by core collapse supernovae [HI]. Those previous multidimensional 
models for core collapse supernovae that did produce explosions tended to greatly exceed 
these limits (see, e.g., [TTJ [HE1 [HE]). To compensate, modelers have been forced to 
posit the fallback of a considerable amount of matter onto the neutron star, occurring 
on a timescale longer than was simulated. While the decreased neutronization seen 
in Boltzmann models reduces the need to invoke fallback, it also makes any fallback 
scenario more complicated, since the most neutron-rich material may no longer be the 
innermost. Because of the impact of the neutrinos on the nucleosynthesis, we fully expect 
the nucleosynthesis products from Boltzmann explosion simulations to be qualitatively 
different, both in composition and spatial distribution, from either parameterized bomb 
[HE] or piston [21] nucleosynthesis models or the present generation of models of the core 
collapse mechanism. 

Figure E] demonstrates the effect of the neutrino interactions on the innermost 
ejecta. In this example, the neutrino interactions of a model similar to [21] have been 
modified in order to produce an explosion [67J, but it is expected that realistic neutrino- 
driven explosions would be similar. In the innermost ejecta, the shock fully dissociates 
the matter, so neutrino interactions with free nucleons dominate, producing a marked 
increase in the electron fraction. In more distant regions, cooler peak temperatures 
will cause more poorly known v and e 1 * 1 interactions with heavy nuclei to dominate. 
These interactions, as well as neutral current inelastic neutrino scattering off these 
nuclei [H2]; are also important to the thermal balance, affecting the a-richness of the 
ejecta and, thereby, the abundance of important nuclei like 44 Ti, 57 Fe, 58 Ni and 60 Zn 
[21] . Thus, there is a clear need for improved neutrino nucleus interaction rates in order 
to accurately calculate the iron-peak nucleosynthesis from core collapse supernovae. 
Because the degree of neutronization is much less than in deeper layers of the star, 
several nuclei of interest are directly accessible via the KARMEN technique: 40 Ca, 45 Sc, 
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Figure 6. The temporal history of the electron fraction, temperature and density 
for a zone in the innermost ejecta of a supernova model. Note that the effects of 
neutrino capture dominate the effects of neutrino emission (e^ capture), decreasing 
the ncutronization. 



51 V, 55 Mn, 59 Co. However, theoretical calculations will still be necessary to provide full 
coverage of the many species present in significant concentrations. 

5.2. Neutrino Nucleosynthesis 

Neutrino nucleosynthesis is driven by the spallation of protons, neutrons, and alpha 
particles from nuclei in the overlying stellar layers by the intense neutrino flux that is 
emanating from the central proto- neutron star powering the supernova [HE]- Moreover, 
neutrino nucleosynthesis continues after the initial inelastic scattering reactions and 
the formation of their spallation products. The neutrons, protons, and alpha particles 
released continue the nucleosynthesis through further reactions with other abundant 
nuclei in the high-temperature supernova environment, generating new rare species. 
Neutrino nucleosynthesis occurs in two stages: (1) through the neutrino irradiation and 
nuclear reactions prior to shock arrival and (2) through the continuation of nuclear 
reactions induced by neutrinos as the stellar layers expand and cool. The suggestion 
has been made jHH] that neutrino nucleosynthesis is responsible for the production of 
11 B, 19 F, and two of Nature's rarest isotopes: 138 La and 180 Ta . 
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Observations of the abundance of boron vary linearly with metallicity, implying 
that primary mechanisms that operate early in the history of our galaxy produce as 
much of these isotopes as secondary (quadratic) mechanisms that operate after the 
Galaxy has been enriched with metals. The competing mechanism for boron formation, 
cosmic ray spallation, is a secondary process. According to current models, neutrino 
nucleosynthesis in supernovae, which is a primary process, favors the production of 11 B 
over 10 B. However laboratory calibration of the spallation channels producing these two 
isotopes is needed. Used in conjunction with future HST observations discriminating 
between 10 B and n B, this measurement would be invaluable in resolving this controversy 
and supporting (or refuting) the suggestion that neutrino nucleosynthesis in supernovae 
is an important source of n B in the Galaxy [HH]. n B and 10 B are produced through the 
following spallation channels: 

12 C(z/, z/p) n B 

12 C(z/, z/'n) 11 C(e + i/) 11 B 

12 C(z/, z/pn) 10 B. 

The final abundance of 19 F produced in a supernova can serve as a "supernova 
thermometer." If the abundance of 19 F produced in the supernova is attributed to 
neutrino nucleosynthesis, the ratio of [ 19 F/ 20 Ne]/[ 19 F/ 20 Ne] Q (the denominator is the 
measured ratio in the Sun) is a measure of the mu and tau neutrinosphere temperatures 
|70j . 19 F is produced through the following spallation channels: 

20 Ne(z/, i/n) 19 Ne(e+i/ c ) 19 F 
20 Ne(z/, u'p) 19 F. 

Recent models [7T], using improved neutrino nucleus reaction rates, show marked 
decreases in the production of 19 F. 

No obvious astrophysical site for the production of the rare isotopes 138 La and 
180 Ta has been proposed. That they can be produced via neutrino nucleosynthesis in 
supernovae is compelling, and may be very important in that their existence, however 
rare, may be a fingerprint of the neutrino process. If so, they potentially provide powerful 
diagnostics of the physics of the outer layers of the supernova. 138 La and 180 Ta are 
produced through the following spallation channels: 

139 La(z/, !/n) 138 La 
181 Ta(^, z/n) 180 Ta. 

Recent models 71 imply significantly larger production of these isotopes, enhancing 
the possibility that these isotopes originate in supernovae. 

Experiments to measure the cross sections for all of these spallation channels are 
worthy of consideration for experimental determinations at stopped pion facilities. 
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5.3. Nucleosynthesis in the Neutrino-Driven Wind 

The astrophysical r-process (rapid neutron capture process) is responsible for roughly 
half of the Solar System's supply of elements heavier than iron. While the nuclear 
conditions necessary to produce the r-process are well established (see, e.g., |Z21), the 
astrophysical site remains uncertain. The leading candidate is the neutrino-driven wind 
emanating from the proto- neutron star after a core collapse supernova is initiated |73j . 
However, other plausible sites have been suggested [31 EH]- As the neutrino-driven 
wind expands rapidly and cools, the nuclear composition is dominated by a-particles 
and free neutrons with a small concentration of iron group nuclei. As temperatures 
continue to drop, charged particle reactions "freeze out" while neutron capture reactions 
continue on the "seed" heavy nuclei present at freeze-out. Neutron capture (n, 7) 
reactions are balanced by their inverse photodisintegration (7, n) reactions, establishing 
an equilibrium between the free neutrons and the nuclei in the wind. Because of the 
high concentration of free nucleons, this (n, 7)-(7, n) equilibrium among isotopes of the 
same element produces nuclei that are quite neutron rich. (3 decays of nuclei with half 
lives that are short compared to the time scale for the r-process link these (n, 7)-(7, n) 
clusters, producing nuclei with higher Z and leading to the synthesis of heavier elements 

Eg. 

Qian et al have demonstrated that neutrino-induced reactions can significantly 
alter the r-process path and its yields in both the (n, 7) <-> (7, n) equilibrium phase and 
the "postprocessing phase" that occurs once these reactions fall out of equilibrium. 
In the presence of a strong neutrino flux, i/ e -mduced charged current reactions on the 
waiting point nuclei at the magic neutron numbers N = 50, 82, 126 might compete 
with beta decays and speed up passage through these bottlenecks. Also, neutrinos can 
inelastically scatter on r-process nuclei via z/ e -induced charged-current reactions and v- 
induced neutral-current reactions, leaving the nuclei in excited states that subsequently 
decay via the emission of one or more neutrons. This processing may for example 
shift the abundance peak at A = 195 to smaller mass. Extending this, Haxton et 
al jZH| pointed out that neutrino postprocessing effects would provide a fingerprint 
of a supernova r-process. Eight abundances are particularly sensitive to the neutrino 
postprocessing: 124 Sn, 125 Te, 126 Te, 183 W, 184 W, 185 Re, 186 W, and 187 Re. Observed 
abundances of these elements are consistent with the postprocessing of an r-process 
abundance pattern in a neutrino fluence consistent with current supernova models. 
On a more pessimistic note, Meyer, McLaughlin, and Fuller [7D J have investigated 
the impact of neutrino-nucleus interactions on the r-process yields and have discovered 
that electron neutrino capture on free neutrons and heavy nuclei (in the presence of a 
strong enough neutrino flux) can actually hinder the r-process by driving the neutrino- 
driven wind proton rich, posing a severe challenge to theoretical models. This push to 
lower neutronization makes the early phases of the neutrino-driven wind a candidate for 
production of the light p-process nuclei like 74 Se, 78 Kr, 84 Sr and 92 Mo jHS]- 

During the r-process and subsequent postprocessing in the supernova neutrino 
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fluence, neutrinos interact with extremely neutron-rich, radioactive nuclei. Thus, 
relevant direct neutrino-nucleus measurements cannot be made. However, indirect 
measurements of charged- and neutral-current neutrino-nucleus interactions on heavy 
stable nuclei for A > 80 would be invaluable as a gauge of the accuracy of theoretical 
predictions. 

6. Supernova Neutrino Detection 

The nineteen neutrino events detected by 1MB and Kamiokande from SN1987A 
confirmed the basic supernova paradigm — that core collapse supernovae mark the 
formation of a neutron star and release copious amounts of neutrinos — and signaled the 
birth of extra-Solar-System neutrino astronomy. For a Galactic supernova, thousands 
of events will be seen by Super-K and SNO, which for the first time will give 
us detailed neutrino "light curves" and bring us volumes of information about the 
deepest regions in the explosion. In turn, these light curves can be used to test 
and improve supernova models, thereby improving predictions about the explosion and 
resultant nucleosynthesis. Moreover, from these detailed neutrino light curves and an 
understanding of the effects of neutrino oscillations, interesting insight could be gained 
about the density structure of the supernova progenitor. Among the neutrino-nucleus 
interactions of relevance for supernova neutrino detection are neutrino interactions on 
2 H, 16 0, 56 Fe and 206,207,208p b _ 

6.1. Deuterium 

Neutrino experiments that use heavy water, like the Sudbury Neutrino Observatory 
(SNO), can detect supernova neutrinos via four main channels: 

e~(v, u')e~ 
d(z/, z/n)p 
d(i/ e ,e"p)p 
&{y C) e + n)n 

Measurement of the reaction d(z/ e ,e~p)p, which is being considered to calibrate the 
reaction p(p,e + i/ e )d (part of the pp chain of reactions powering the Sun), would also 
provide a calibration for heavy water neutrino detectors. Monte Carlo studies suggest 
that for the source brightness predicted for the ORNL SNS, two years of data in 
approximately thirty fiducial tons of D2O would yield a cross section measurement 
with an accuracy of a few percent [HO] , which in turn will enable a more accurate 
interpretation of the SNO data from the next Galactic supernova. 

6. 2. Oxygen 

The charged-current reaction 16 0(z/ e , e~) 16 F is the principle channel for electron neutrino 
interactions for thermal sources in the range T Uc > 4 — 5 MeV, and its rate exceeds 



Supernova Science at Spallation Neutron Sources 



19 



that of neutrino-electron scattering by an order of magnitude for T Vc > 7 — 9 MeV 
[HT] . Moreover, the electron angular distribution is strongly correlated with the 
electron neutrino energy, providing a way to measure the incident neutrino energy and, 
consequently, the electron neutrino spectra. By inference, one would then be able to 
measure, for example, the electron neutrinosphere temperature, providing an additional 
supernova thermometer [5U] . 

In addition, the appearance of back-angle electron emission from this reaction in, for 
example, Super-K would result from energetic electron neutrinos, more energetic than 
predicted by supernova models, providing further evidence for flavor oscillations and 
thereby information about the mu and tau neutrino spectra emanating from supernovae 
[5Uj . Mu and tau neutrinos in the stellar core couple to the core material only via 
neutral currents, whereas electron neutrinos and antineutrinos couple via both neutral 
and charged currents. As a result, the former decouple at higher density and, therefore, 
temperature, and have harder spectra. Utilizing reactions on 16 0, Langanke, Vogel, 
and Kolbe [H2] have suggested a novel way of also unambiguously identifying mu and 
tau neutrino signatures in Super-K. The larger average energies for these neutrino 
flavors may be sufficient to excite giant resonances via the neutral-current reactions 
16 0(V^ r , z/^ T ) 16 0*. These resonances are above particle threshold and subsequently 
decay via the emission of protons, neutrons, and gamma rays. The gamma rays would 
provide the mu and tau neutrino signatures. The two decay channels are: 16 0*(, 7n) 15 
and 16 0*(, 7p) 15 N. However, potential channels for observing the mu and tau neutrinos 
from supernovae must be reexamined in light of recent work (see, e.g., IM| 
which indicates that nucleon-nucleon bremsstrahlung and the effects of nuclear recoil 
in neutrino-nucleon scattering significantly soften the mu and tau neutrino spectra, 
lessening their energy excess over electron neutrinos (see Figure |7J). 

Thus, accurate measurements of both charged- and neutral-current neutrino cross 
sections on 16 would serve as a foundation for interpreting the neutrino data from the 
next Galactic core collapse supernova and for using the data to potentially observe the 
mu and tau neutrino spectra as it is emitted from the proto-neutron star. An experiment 
to measure the cross section for: 

16 0(z/ c ,e-) 16 F 

should be a high priority for a stopped pion facility. Further useful experiments could 
focus on the cross sections for: 

16 0(^,z/n 7 ) 15 
16 0(^,< P7 ) 15 N. 

6.3. Iron and Lead 

The use of iron and lead in supernova neutrino detectors like the proposed OMNIS 
detector would provide another way of detecting the mu and tau neutrino spectra in core 
collapse supernovae |89j . Iron has a sufficiently high threshold for neutron production via 
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Figure 7. Comparison of the neutrino spectra from supernova simulations using 
the standard ^B] opacities (left) and updated opacities (right; including the effects 
of neutino-nucleon absorption and elastic scattering 86 , neutrino-nucleon inelastic 
scattering and bremsstrahlung 87 , and weak magnetism [EE])- The simulations, 
initiated from a 13 M progenitor, are fully general relativistic, and the spectra are 
computed at a radius of 500 km, 100 milliseconds after bounce |85j . 



charged- current neutrino interactions that such production is negligible, whereas in lead 
neutrons are produced by both charged- and neutral-current interactions. Oscillations 
between the more energetic mu and tau neutrinos and the electron neutrinos would boost 
the charged- current event rate while leaving the neutral-current rate roughly unchanged. 
Thus, the ratio of the event rate in lead to that in iron would serve as a further constraint 
on the extent of neutrino oscillations and the emitted mu/tau spectra. However, this 
potential channel for observing the mu and tau neutrinos must also be reexamined in 
light of the softening of the mu and tau neutrino spectra. 

To further the development of a detector like OMNIS, experiments to measure 
the neutrino-iron and neutrino-lead cross sections have been proposed. For iron, the 
neutral-current reaction: 

56 Fe(z/, z/n) 55 Fe 

dominates. For lead, a total cross section would be measured resulting from the following 
neutral- and charged-current channels: 

208 Pb(u, z/n) 207 Pb 
208 Pb(z/, z/'2n) 206 Pb 
208 Pb( I / e ,e-n) 207 Bi 

including also the channels for the isotopes 206 Pb and 207 Pb. For this reason, as well as 
their importance to nucleosynthesis, iron and lead cross section measurements should 
be among the first experiments at a stopped pion facility. 
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Figure 8. The expected neutrino spectra from the Spallation Neutron Source at 
ORNL. 



7. Conclusion 



A new facility to measure neutrino-nucleus cross sections, such as those feasible 
at stopped pion facilities, would provide an experimental foundation for the many 
neutrino-nucleus weak interaction rates needed in supernova models. With a neutrino 
source as intense as the ORNL SNS, for example, we are presented with a unique 
opportunity, given the overlap between the facility and supernova neutrino spectra 
(compare Figure [7| & EJ) , to make such measurements. This would enable more realistic 
supernova models and allow us to cull fundamental physics from these models with 
greater confidence when their predictions are compared with detailed observations. 
Charged- and neutral-current neutrino interactions on nuclei in the stellar core play 
a central role in supernova dynamics and nucleosynthesis and are also important for 
supernova neutrino detection. Measurements of these reactions on select, judiciously 
chosen targets would provide an invaluable test of the complex theoretical models used 
to compute the neutrino-nucleus cross sections. 



Supernova Science at Spallation Neutron Sources 



22 



Acknowledgments 

The work has been partly supported by NASA under contract NAG5-8405, by the 
National Science Foundation under contract AST-9877130, by the Department of 
Energy, through the PECASE and Scientic Discovery through Advanced Computing 
Programs, and by funds from the Joint Institute for Heavy Ion Research. A.M. is 
supported by the Oak Ridge National Laboratory, which is managed by UT-Battelle, 
LLC, for the U.S. Department of Energy under contract DE-AC05-00OR22725. The 
authors would like to acknowledge many illuminating discussions with Frank Avignone, 
John Beacom, Jeff Blackmon, Dick Boyd, David Dean, Yuri Efremenko, Jon Engel, 
George Fuller, Wick Haxton, Ken Lande, Karlheinz Langanke, Matthias Liebendorfer, 
Gabriel Martinez-Pinedo, Gail McLaughlin, Mike Strayer, and Friedel Thielemann, all 
of which contributed significantly to this manuscript. 

References 

[1] Burrows A and Lattimcr J M 1985 ApJ 299 L19-L22 

[2] Wilson J R 1985 Numerical Astrophysics ed J M Centrella J M LcBlanc and R L Bowers (Boston: 

Jones and Bartlett) pp 422-434 
[3] Bcthc H H and Wilson J R 1985 ApJ 295 14-23 

[4] Langanke K and Martinez-Pinedo G 2000 Nucl. Phys. A 673 481-508 
[5] Fuller G M, Fowler W A, and Newman M J 1985 ApJ 293 1-16 
[6] Baron E, Cooperstein J, and Kahana S 1985 Phys Rev. Lett. 55 126-129 
[7] Lattimer J M and Swesty F D 1991 Nucl. Phys. A 535 331-376 

[8] Aufderheide M B, Fushiki I, Woosley S E, and Hartmann D H 1994 ApJS 91 389-417 

[9] Fryer C, Burrows A, and Benz W 1998 ApJ 496 333-351 
[10] Wang L, Howell D A, Hoflich P, and Wheeler J C 2001 ApJ 550 1030-1035 
[11] Janka H-T and Miiller E 1996 A&A 306 167-98 
[12] Burrows A and Goshy J 1993 ApJ 416 L75-L78 
[13] van Riper K A and Lattimcr J M 1981 ApJ 249 270-289 
[14] Cooperstein J, van den Horn L J, and Baron E A 1986 ApJ 309 653-666 
[15] Arnett W D 1977 ApJ 218 815-833 
[16] Bruenn S W 1985 ApJS 58 771-841 

[17] Myra E S, Bludman S A, Hoffman Y, Lichenstadt I, Sack N, and van Riper K A 1987 ApJ 318 
744-759 

[18] Wilson J R and Mayle R W 1993 Phys. Rep. 227 97-111 
[19] Janka H-T 1992 A&A 256 452-458 

[20] Messer O E B, Mezzacappa A, Bruenn S W, and Guidry M W 1998 ApJ 507 353-60 
[21] Mezzacappa A, Liebendorfer M, Messer O E B, Hix W R, Thielemann F-K, and Bruenn S W 2001 
Phys. Rev. Lett. 86 1935 

[22] Liebendorfer M, Mezzacappa A, Thielemann F-K, Messer O E B, Hix W R, and Bruenn S W 2001 

Phys. Rev. D 63 103004 
[23] Nomoto K and Hashimoto M 1988 Phys. Rep. 163 13 
[24] Woosley S E and Weaver T A 1995 ApJS 101 181 

[25] Heger A, Langanke K, Martinez-Pinedo G, and Woosley S E 2001 Phys. Rev. Lett. 86 1678-1681 
[26] Rampp M and Janka H-T 2000 ApJ 539 L33-L36 

[27] Liebendorfer M, Mezzacappa A, and Thielemann F 2001 Phys. Rev. D 63 19643 
[28] Mezzacappa A and Messer O E B 1999 J. Comp. Appl. Math 109 281-319 



Supernova Science at Spallation Neutron Sources 



23 



Liebendorfer M, Messer O E B, Mezzacappa A, and Hix W R 2001 Proceedings of the 20th Texas 
Symposium on Relativistic Astrophysics ed J C Wheeler and H Martel (Melville: AIP) pp 472- 
477 

Oda T, Hino M, Muto K, Takahara M, and Sato K 1994 Atomic Data and Nuclear Data Tables 
56 231 

Mazurek T J 1973 Ph.D. thesis Ycshiva University 
Hansen C J 1966 Ph.D. thesis Yale University 
Yahil A 1983 ApJ 265 1047-1055 

Messer O E B, Liebendorfer M, Mezzacappa A, Hix W R, and Bruenn S W 2003 ApJ in preparation 
Mezzacappa A and Bruenn S W 1993 ApJ 405 637-668 
Liebendorfer M, Rosswog S, and Thielemann F-K 2002 ApJS 141 229-246 

Liebendorfer M, Messer O E B, Mezzacappa A, Bruenn S W, Cardall C Y, and Thielemann F-K 

2003 ApJS in press 
Mezzacappa A and Bruenn S W 1993 ApJ 405 669-684 
Fuller G M 1982 ApJ 252 741-764 

Cooperstcin J and Wambach J 1984 Nucl. Phys. A 420 591-620 
Langanke K, Kolbc E, and Dean D J 2001 Phys. Rev. C 63 32801 

Langankc K, Martinez-Pinedo G, Sampaio J M, Dean D J, Hix W R, Messer O E B, Mezzacappa 

A, Liebendoerfer M, Janka H-T, and Rampp M 2003 Phys. Rev. Lett. 90 241102 
Hix W R, Messer O E B, Mezzacappa A, Liebendoerfer M, Sampaio J M, Langanke K, Martinez- 
Pinedo G, and Dean D J 2003 Phys. Rev. Lett, in press 
Langanke K, Martmcz-Pincdo G, and Sampaio J M 2001 Phys. Rev. C 64 55801 
Clifford F E and Taylcr R F 1965 MNRAS 129 104 
Hartmann D, Woosley S E, and El Eid M F 1985 ApJ 297 837-845 
Hix W R 1995 Ph.D. thesis Harvard University 
Bravo E and Garcfa-Scnz D 1999 MNRAS 307 984-992 
El Eid M F and Hillebrandt W 1980 A&AS 42 215-226 

Brachwitz F, Dean D, Hix W, Iwamoto K, Langanke K, Martinez-Pinedo G, Nomoto K, Strayer 

M R, and Thielemann F-K 2000 ApJ 536 934-947 
Maschuw R 1998 Prog. Part. Phys. 40 183 

Mezzacappa A, Calder A C, Bruenn S W, Blondin J M, Guidry M W, Strayer M R, and Umar 

A S 1998 ApJ 493 848-62 
Keil W, Janka H-T, and Miiller E 1996 ApJ 473 L111-L114 

Buras R, Rampp M, Janka H-T, and Kifonidis K 2003 Phys. Rev. Lett. 90 241101 
Bruenn S W and Dincva T 1996 ApJ 458 L71-L74 

Hcrant M, Bcnz W, Hix W R, Fryer C L, and Colgate S A 1994 ApJ 435 339-361 
Burrows A, Hayes J, and Fryxcll B A 1995 ApJ 450 830-850 

Swesty F D 1998 Stellar Evolution, Stellar Explosions, and Galactic Chemical Evolution ed A 

Mezzacappa (BristoklOP) p 539-48 
Mezzacappa A, Calder A C, Bruenn S W, Blondin J M, Guidry M W, Strayer M R, and Umar 

A S 1998 ApJ 495 911-26 
Fryer C L and Warren M S 2002 ApJ 574 L65-LL68 
Blondin J, Mezzacappa A, and DeMarino C 2003 ApJ 584 971-980 
Bruenn S W and Haxton W C 1991 ApJ 376 678-700 

Hoffman R D, Woosley S E, Fuller G M, and Meyer B S 1996 ApJ 460 478-88 
Trimble V 1991 A&AR 3 1-46 

Kifonidis K, Plewa T, Janka H-T, and Miiller E 2000 ApJ 531 L123-L126 
Thielemann F-K, Nomoto K, and Hashimoto M 1996 ApJ 460 408-36 

Martinez-Pinedo G, Hauser P, Hix W, Liebendrfer M, Mezzacappa A, and Thielemann F-K 2003 
ApJ in preparation 

[68] Woosley S E, Hartmann D H, Hoffman R D, and Haxton W C 1990 ApJ 356 272-301 



Supernova Science at Spallation Neutron Sources 



24 



Hartmann D, Myers J, Woosley S, Hoffman R, and Haxton W 1999 ASP Conf. Scr. 171: LiBeB 

Cosmic Rays, and Related X- and Gamma-Rays p 235 
Haxton W 1999 in Neutrinos in Physics and Astrophysics From 1CT 33 to 10 28 cm (TASI 98) ed P 

Langacker (Singapore:World Scientific) p 432-87 
Heger A, Kolbe E, Haxton W, Langanke K, Martinez-Pinedo G, and Woosley S 2003 Phys. Rev. 

Lett, submitted 

Kratz K-L, Bitouzct J-P, Thielemann F-K, Moller P, and Pfciffcr B 1993 ApJ 403 216-238 
Woosley S E, Wilson J R, Mathews G J, Hoffman R D, and Meyer B S 1994 ApJ 433 229-246 
Wheeler J C, Cowan J J, and Hillebrandt W 1998 ApJ 493 L101-L104 
Freiburghaus C, Rosswog S, and Thielemann F-K 1999 ApJ 525 L121-L124 
Meyer B S 1994 ARA&A 32 153-190 

Qian Y-Z, Haxton W C, Langanke K, and Vogel P 1997 Phys. Rev. C 55 1532-1544 

Haxton W C, Langanke K, Qian Y-Z, and Vogel P 1997 Phys. Rev. Lett. 78 2694-2697 

Meyer B S, McLaughlin G C, and Fuller G M 1998 Phys. Rev. C 58 3696-3710 

Workshop on neutrino-nucleus physics using a stopped-pion facility (Orland) 2000 Tech report 

Oak Ridge National Laboratory 
Haxton W C 1987 Phys. Rev. D 36 2283-2292 

Langanke K, Vogel P, and Kolbe E 1996 Phys. Rev. Lett. 76 2629-2632 
Thompson T A, Burrows A, and Horvath J E 2000 Phys. Rev. C 62 035802 
Raffelt G G 2001 ApJ 561 890-914 
Bruenn S W 2004 ApJ in preparation 

Reddy S, Prakash M, and Lattimcr J M 1998 Phys. Rev. D 58 013009 
Hannestad S and Raffelt G 1998 ApJ 507 339-352 
Horowitz C J 2002 Phys. Rev. D 65 43001 

Boyd R N, Murphy A S, and OMNIS collaboration 2001 Nucl. Phys. A 688 386-389 



